Quasi-degenerate self-trapping in One-Dimensional Charge Transfer Exciton 
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The self-trapping by the nondiagonal particle-phonon interaction between two quasi-degenerate 
energy levels of excitonic system, is studied. We propose this is realized in charge transfer exciton, 
where the directions of the polarization give the quasi-degeneracy. It is shown that this mechanism, 
unlike the conventional diagonal one, allows a coexistence and resonance of the free and self-trapped 
states even in one-dimensional systems and a quantitative theory for the optical properties (light 
absorption and time-resolved luminescence) of the resonating states is presented. This theory gives 
a consistent resolution for the long-standing puzzles in quasi-one-dimensional compound A-PMDA. 
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The charge transfer (CT) exciton in one-dimension 
(ID) attracts great interests due to its pecuhar features, 
e.g., strong interaction with the phonons, magnons, and 
large nonlinearity in the optical responses. However the 
theoretical understanding of it is rather limited compared 
with that of Wannier and Frenkel excitons, because of its 
intermediate radius and strong coupling nature. E.g., al- 
though the prototypical compound Anthracene-PMDA 
(A-PMDA), which shows unusual vibronic structure in 
the optical spectra, is thoroughly experimentally inves- 
tigated from 70-th till nowdays even qualitative 
theoretical understanding it's optical response is missing. 

The most dramatic phenomenon which occurs at in- 
termediate or strong coupling with the phonons is the 
self-trapping where the free (F) state of polaron with 
small lattice deformation coexists and resonates with the 
state with strong lattice relaxation PpO) . Although the 
latter state is Bloch invariant it is traditionally classi- 
fied as being self-trapped (ST) to stress larger degree of 
ionic distortion around the quasiparticle. However, this 
self-trapping scenario was not considered yet since the F- 
ST resonance is prohibited by Rashba-Toyozawa theorem 
(RTT) in ID systems Ahhough the RTT is proba- 

bly valid for one quasiparticle state in the relevant energy 
range, the neglect of it's domain of definition caused se- 
vere delusion in the physics of ID systems because RTT 
was treated as the strict ban for F-ST coexistence in ID. 
However, when there are two quasi-degenerate quasipar- 
ticle states, i.e. the conditions of the theorem are vio- 
lated, one can consider the quasi-degenerate self-trapping 
(QDST) mechanism, when the F-ST resonance is driven 
by nondiagonal interaction with phonons with respect to 
quasiparticle levels. An example of QDST was observed 
in the mixed valence systems [ |T]] |. 

In this paper we first demonstrate by exact diagram- 
matic Monte Carlo method that QDST is general 
phenomenon pertinent for ID lattice, study the optic 
properties of resonating QDST states, and resolve long- 
standing puzzle of the quasi-lD system A-PMDA K 



QDST existence in ID by exact solution. The quasi-lD 
compound A-PMDA consists of one-dimensional arrays 
of alternating donor (D) and acceptor (A) molecules [|l|,D. 
There are two possible configurations of CT exciton, i.e., 
A'^D^A^ and A^D'^A'^, and the symmetric (j)s{r) and 
antisymmetric (jfa (r) linear combination of these two con- 
stitutes the nearly degenerate two branches of excitons. 

The minimal model to demonstrate F-ST coexistence 
in ID systems involves one optic phonon branch with fre- 
quency Lo — 0.1, two quasiparticle branches with energies 
£1,2(9) = Ai^2 + 2[1 — cos{q)], where Ai = and A2 — 1, 
and the quasiparticle-phonon coupling 
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Here Cj,k and bq are annihilation operators for the quasi- 
particle of branch j with momentum k and for the 
phonon with momentum q. One can restrict himself 
by one phonon mode since, even in case of the differ- 
ent symmetry of excitons, the same phonon branch con- 
tributes both to the diagonal and off-diagonal coupling 
|l3| and drop, without affecting the conclusions, the q- 
dependence {4>ij{q) = 1) of the interaction for simplic- 
ity. The energy E^, the average number of phonons 
N = 'l2q^q^q I ^^^d Z-factors of the n-phonon 
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(Here | vac) is the vacuum state.) were calculated by the 
exact diagrammatic Monte Carlo technique which 
was generalized in the present paper to the case of two- 
quasiparticle bands with interaction Hamiltonian (fil) . In 
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FIG. 1. Energy (solid line, circles) and iV (dashed line, 
squares) dependence on 712 at 711 — and 711 = 0.05. Inset: 
Z-factors distribution for 712 — 0.2087. The statistic errorbars 
are much smaller than the point size. 



the vicinity of critical 7^3 ~ 0.2087 the dependence of the 
energy Ej^ on 712 rapidly changes, and the average num- 
ber of phonons N demonstrates sharp crossover to higher 
values (Fig. |l|). Hence, according to HJl^, the F-state, 
which is lowest for small values of nondiagonal interaction 
constant, resonates with the ST-state (see inset in Fig. |^) 
at 7^2 and ST-state is the lowest one for large values of 
712 10 • The system indeed undergoes crossover around 
the critical constant rather than a phase transition due 
to mixing of resonating eigenstates . 

Strong- coupling limit model. If the quasiparticle is ex- 
citon, one can study the optical transition between the 
ground state G and the exciton-polaron states with 
zero momentum (see Fig. To study the optical prop- 
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FIG. 2. Schematic configurational diagram. The dotted ar- 
rows correspond to nonradiative processes, and the dashed ar- 
rows indicate the radiative (absorption/emission) transitions. 

erties we consider the case when the coexistence occurs 



in the strong coupling limit, i.e. when 7x2/'^ ^ ^ and, 
therefore, the adiabatic approximation is justified [ p^ . 
In this case one can find for given lattice deformation 
Q = {Qi, ---jQn} (N is the number of lattice modes) 
the adiabatic potential relieves eA(Q) (A is the index of 
lower I and upper u sheet) and express the eigenfunc- 
tions il}\{r,Q) = a^((5)0i(r) -I- (Q)?!'2(i") in terms of 
the excitonic wave functions 0i,2(r) for undistorted lat- 
tice . In case of F-ST coexistence the lower sheet 
of the adiabatic potential contains two minima at dif- 
ferent lattice distortions Q* {t = ST) which corre- 
spond to states with different lattice relaxation energy. 
When the minima are separated by large potential barrier 
one can neglect the Q-dependence of the electronic func- 
tions and introduce vibrational states xV^iQ ~ Q*) = 
n^Lixf"(Qn ~ Qn) with quantum numbers {/„}. The 
low laying eigenstates of the adiabatic approximation are 



nxP(Q-Q*); t 



F,ST, (2) 



with the eigenvalues E^^ = e;('9*) + E,T=i ^r(/n + l/2) 
The adiabatic basis (|^) is an approximate one due to non- 
diagonal matrix elements of nonadiabatic operator [ pT| . 
E.g., at critical 7^2 the ground state is already a mix- 
ture of weakly/strongly deformed states (inset in Fig. |^). 
However, the basis (^ is valid when interaction constant 
is close enough to the critical value 7^2 to create two wells 
but far enough to neglect the mixing of F and ST sets. 

Spectrum of optical absorption. The optical proper- 
ties of A-PMDA give the evidence of the strong exciton- 
phonon interaction The low energy part of the ab- 
sorption spectrum consists of the zero phonon line (ZPL) 
and a broad band with the full width at half maxi- 
mum (FWHM) of around |,| 500 cm . Therefore, 
the optical response qualitatively remind lineshape of 
F-color centers absorption spectra which correspond to 
transitions between displaced oscillators How- 
ever, all attempts to interpret the optical properties of 
A-PMDA [0_0 encountered apparent contradictions with 
the model | |l7| , [l8t : (i) absorption and emission spec- 
tra can not be described as one series of the equidis- 
tant lines (ii) the value of Huang- Rhys parameter 
S, evaluated from the ratio of ZPL and total oscilla- 
tor strength, is inconsistent with the intensity distribu- 
tion of the broad band [E.g., the value 5* = 4 extracted 
from the oscillator strengths of ZPL leads to the estimate 
FWHM=140cm~i. |].]. 

These long-standing puzzles can be interpreted (such 
treatment has never been attempted due to the RTT) 
as response of the coexisting F and ST states. To cal- 
culate induced cross section a'^^^{E) at energy E of op- 
tical absorption which corresponds to transitions from 

the set = Mr)xh"HQ) of xh^HQ) phonon 

states of ground electronic state 0G(r) to the set of ex- 
cited states in one has to average over initial states 
{a} and sum over final states in (0). Using the standard 
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approach defining the electronic dipole matrix ele- 
ments Mi = J dr(0G(r))*r</)j(r), and introducing an av- 
erage optic phonon frequency Ut one can take advantage 
of Huang-Rhys method [ pT| and obtain for zero temper- 
ature 



E 

=F,ST 



MtF^^E), 



(3) 



where 



-St 



Sf 



(J 7+ ' 7)! 



p=0 



(Here 5* = {2N)~^u,t J2n=iiQi? is the Huang-Rhys fac- 
tor.) is the normalized shape function, and 



Mt = 



E 



(M,aKQ*))*M,4(g*; 



(5) 



is the weighting coefficient. Here e and c is electron 
charge and light velocity, (E) is the average energy of 
transitions, TP) = (Tf /2^)(x2-|-(T^/2)2)-i, and Tf 
is the spectroscopic linewidth. We note that for T ^ 
one can also use approach since the average is per- 
formed over the one well of the ground state pO[ |. 

The contradiction (i) in the optical absorption spec- 
trum of A-PMDA is resolved when one note that experi- 
mental data 1^ contain two series of satellites in absorp- 
tion spectrum: one arises from the ZPL with the step 
LUF ^27 cm-i (0 - 29 - 55 - 83 - 111 cm'^) and the 
second one (13 - 29 - 47 - 62 - 77 - 88 cm^^) with the 
overtone u)st ^ 12 cm~^ originates from the 13 cm^^ 
satellite. The contradiction (ii) gains natural explana- 
tion since the spectrum is the superposition of the two 
contributions (^ from F and ST states with different pa- 
rameters Sf and Sst (Fig. ||). 
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FIG. 3. Comparison of light absorption experimental data 
(diamonds) with model curve (solid line) consisting of optic 
response of coexisting F (dashed line) and ST (dotted line) 
states with Sf = 2.5, Sst = 13.3, and Mf/Mst = 0.82. 
The onsets are eF,ST(Q^'®^) = 18322/18334 cm"^ and 
linewidthes are T f,st = 29/72 cm~^ (T^T" = 2.5 cm"^). 



Luminescence from QDST states. The main puzzle in 
luminescence spectra of A-PMDA is that although the 
pattern of satellites is symmetrically repeated in absorp- 
tion and luminescence spectra, the intensity distribution 
of phonon sidebands is significantly different ||l^ . 

The experimental conditions of time-resolved lumines- 
cence spectroscopy are equivalent to the situation when 
at a moment i = the system is excited to the up- 
per sheet U whereas F and ST wells are empty. The 
rates of the intrawell nonradiative processes are signifi- 
cantly larger than that of interwell ones wjj^st, wjj^f 
and wsT^F (see Fig. ^), because successive relaxation 
of the system within one well does not demand many- 
phonon transitions whereas those are compulsory for the 
interwell jumps. Due to this first link of the "bottle- 
necks hierarchy" (BH), the states in one well can be 
treated as being in quasi-thermal equilibrium and only 
summary populations Pu{t), Psrit) and PF{t) have to 
be considered as dynamical variables with initial condi- 
tions Puit = 0) = 1 and Psrit = 0) = PF{t = 0) = 0. 
The differential intensities I^'^"^'^{t) of luminescence can 
be expressed in terms of radiative rates Ru, i?ST, and 
Rf as I^'^'^'^{t) = Ru,ST,FPu,ST,F{t), and non time- 
resolved experimental setup is defined by total yields 
^int'^^"^ = Ru,ST.,F Pu,ST,F{t). If the lattice dis- 
tortion in the upper sheet minimum and that for 
ST-minimum of the lower sheet almost coincide whereas 
the lattice relaxation for the minimum F differs 
significantly (such situation can be easily obtained in 
strong-coupling limit model), there is the second link of 
the BH wjj^ST ^ T[iax{wu<^F,wsT^F}, which is due 
to Franck-Condon suppression of nonradiative processes 
U F and ST <-> F. At low temperatures (when, e.g., 
ST state is metastable and the rates wu^st, wu^f, and 
wsT-r-F are negligibly small) the dynamics of the popu- 
lations is reduced to the system of equations |21 



Pu{t) = 
Psrit) = 
Mt) = 



-{WST->F - 
-RFPF{t) 



Wu~.F)Pu{t) 

- RsT)PsT{t) + wu^srPuit) 

- wu^FPu{t) + wsT^FPsrit). 



The differential intensities of the luminescence from 
the states T and F are 

Rstwu^st 



IsAt) 



RST + WST^F ~ WU~ 

-{wu-,ST + 'Wu-,F)t _ 



X 

>ST — Wu^F 

^ — {wsT-^F+RsT)t 



(6) 



and 



lF{t)^RF 



{wij^F + T^) 



[Rf - Wu-,ST - wu^f) 



V 



{Rf — RsT — wst^f) 



(7) 
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respectively. Here 2? = wu^STWsT^F/iRsT + vJST^F — 
wu^ST — wu^f)- Introducing the dimensionless quan- 
tities Vu^F.,U~^ST = Wu^F.U^ST / iwu->ST + Wij^p) 

and rjsT = Rst /wst^f one obtain the integral yields 
as YsT = ?7stVc/^st(1 + Vst)'^ and Yp = (1 + 
VSTVu~*F)i^ + Vst)^^, respectively. There are some 
domains of model parameters when one can observe 
vivid qualitative features. E.g., second link of the BH 
leads to relation Yst/Yf ^ 1 provided the inequahty 
RsT > WST^F takes place. Therefore, due to specific 
BH of the relaxation processes it is possible that, e.g., 
the F state is distinctly observed in the light absorp- 
tion spectra but is almost not seen in the spectra of non 
time-resolved luminescence. At the same time if crite- 
rion Rp < RsT is satisfied, at large times t > Rp^ the 
ratio of differential intensities IsT{t)/lF{t) ^ 1 is oppo- 
site. Note, that this relation implies that A4st > Mf 
because the rate of radiative processes is governed by 
the matrix elements (|^) . Due to the first link of the BH 
at T = only ground vibrational states of the wells ST 
and F are populated at a moment of radiative transition. 
Therefore, the many-well structure leaves intact Huang- 
Rhys approach since no average over the initial states is 
necessary. Hence, the differential spectrum of the lumi- 
nescence is cr(£;, t) = IsT{t)F§^{E)+lFit)F^'^{E). Here 
the emission normalized shape function F™(E) can be 
obtained from by replacing p —p in the function 
S. For the lineshape of non time-resolved experiment one 
has to use integral yields Yst and Yp instead of differ- 
ential intensities (||) and (^. Note, that parameters Sp 
and SsT are the same as in absorption experiment. 

Finally, the difference between absorption and lumines- 
cence spectra pattern arises in A-PMDA due to the BH 
which manifests itself in the mirror image of sidebands 
coming from ST and F series in absorption and lumi- 
nescence spectra, which relative intensities are different 
in two types of experiment. Moreover, due to the BH 
both 13 cm~^ peak Q and the whole broad structure 
|p|,pi)| have inverse relative spectral weights in lumines- 
cence and absorption spectrum. More evidences can be 
found in experiments on time-resolved luminescence be- 
cause proper fit (Fig. ^ of the experimental curve ||] 
demands Aip/Aip < 1 in (||), and, hence, criterion 
Rp < RsT is satisfied. 

Conclusions. We have demonstrated that the coex- 
istence of the F and ST states is possible even in ID 
lattices, owing to the possibility of the QDST mecha- 
nism. When there are more than one electronic levels of 
a particle within the energy range of the lattice relax- 
ation energy, F and ST states can coexist and even res- 
onate. Hence, although the Rashba-Toyozawa theorem 
is not overthrown and remains valid within its domain 
of definition, there is the possibility of the coexisting of 
ST and F states in ID systems. We have found several 
features which are unique for QDST states and can be 
checked in optical experiments. We have shown, that the 



QDST mechanism provides comprehensive description of 
the long-standing mystery of the optical properties of the 
quasi-lD compound A-PMDA. 
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